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Abstract. We study the minimum Manhattan network problem, which is defined as follows. Given a 

set of points called terminals in R d , find a minimum-length network such that each pair of terminals 

'— is connected by a set of axis-parallel line segments whose total length is equal to the pair's Manhattan 

Qh (that is, Li-) distance. The problem is NP-hard in 2D and there is no PTAS for 3D (unless P=AfP). 

•^U Approximation algorithms are known for 2D, but not for 3D. 

We present, for any fixed dimension d and any e > 0, an 0(rf )-approximation algorithm. For 3D, we 
also give a 4(fc — l)-approximation algorithm for the case that the terminals are contained in the union 
of k > 2 parallel planes. 
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1 Introduction 
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In a typical network construction problem, one is given a set of objects to be interconnected such 
CSJ that some constraints regarding the connections are fulfilled. Additionally, the network must be 

of little cost. For example, if the objects are points in Euclidean space and the constraints say 
that, for some fixed t > 1, each pair of points must be connected by a path whose length is 
(^ bounded by t times the Euclidean distance of the points, then the solution is a so-called Euclidean 

t-spanner. Concerning cost, one usually requires that the total length of the network is proportional 
to the length of a Euclidean minimum spanning tree of the points. Such low-cost spanners can be 
constructed efficiently [2]. 

—i In this paper, we are interested in constructing 1-spanners, with respect to the Manhattan 

(or Li-) metric. Rather than requiring that the total length of the network is proportional to the 
minimum spanning tree of the points, our aim is to minimize the total length (or weight) of the 
network. Note that the Euclidean 1-spanner of a set of points is simply the complete graph (if no 
three points are collinear) and hence, its weight is completely determined. Manhattan 1-spanners, 
in contrast, have many degrees of freedom and vastly different weights. 

More formally, given two points p and q in d-dimensional space M rf , a Manhattan path con- 
necting p and q (a p—q M-path, for short) is a sequence of axis-parallel line segments connecting p 
and q whose total length equals the Manhattan distance between p and q. Thus an M-path is 
a monotone rectilinear path. For our purposes, a set of axis-parallel line segments is a network. 
Given a network N, its weight \\N\\ is the sum over the lengths of its line segments. A network N 
Manhattan-connects (or M-connects) two given points p and q if it "contains" a p—q M-path w. 
Note that we slightly abuse the notation here: we mean pointwise containment, that is, we require 
U n — [JN. Given a set T of points — called terminals — in K , a network N is a Manhattan net- 
work (or M-network) for TUN M-connects every pair of terminals in T. The minimum Manhattan 
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(a) an M-network for 
T = {s,s',t,t'} 
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(b) a minimum M- 
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(c) a minimum M- 
network in 3D 



Fig. 1: Examples of M-networks in 2D and 3D. 
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network problem (MMN) consists of finding, for a given set T of terminals, a minimum-weight 
M-network. For examples, see Fig. [TJ 

M-networks have important applications in several areas such as VLSI layout and computational 
biology. For example, Lam et al. |12j used them in gene alignment in order to reduce the size of 
the search space of the Viterbi algorithm for pair hidden Markov models. 



1.1 Previous work 

The 2D-version of the problem, 2D-MMN, was introduced by Gudmundsson et al. [9] . They gave 
an 8- and a 4-approximation algorithm. Later, the approximation ratio was improved to 3 [3 8j and 
then to 2, which is currently the best possible. It was achieved in three different ways: via linear 
programming [5], using the primal-dual scheme [15] and with purely geometric arguments [10]. The 
last two algorithms run in 0{n log n) time, given a set of n points in the plane. A ratio of 1.5 was 
claimed [16], but apparently the proof is incomplete [5]. Chin et al. [6] finally settled the complexity 
of 2D-MMN by proving it NP-hard. 

A little earlier, Muhoz et al. |14j considered 3D-MMN. They showed that the problem is NP-hard 
and that it is NP-hard to approximate beyond a factor of 1.00002. For the special case of 3D-MMN, 
where any cuboid spanned by two terminals contains other terminals or is a rectangle, they gave 
a 2a-approximation algorithm, where a denotes the best approximation ratio for 2D-MMN. They 
posed the design of approximation algorithms for general 3D-MMN as an open problem. 



1.2 Related problems 



As we observe in Section 2.3, MMN is a special case of the directed Steiner forest problem (DSF). 
More precisely, an instance of MMN can be decomposed into a constant number of DSF instances. 
The input of DSF is an edge-weighted directed graph G and a set of vertex pairs. The goal is to 
find a minimum-cost subgraph of G (not necessarily a forest) that connects all given vertex pairs. 
Recently, Feldman et al. [7] reported, for any e > 0, an 0(n 4 ' 5+£ )-approximation algorithm for 
DSF, where n is the number of vertices of the given graph. This bound carries over to dD-MMN. 

An important special case of DSF is the directed Steiner tree problem (DST). Here, the input 
instance specifies an edge- weighted digraph G, a root vertex r, and a subset S of the vertices of G 
to which r must connect. An optimum solution for DST is a minimum- weight r-rooted subtree of G 
spanning S. DST admits an 0(n £ )-approximation for any e > [5J. 

A geometric optimization problem that resembles MMN is the rectilinear Steiner arborescence 
problem (RSA). Given a set of points in R with non-negative coordinates, a rectilinear Steiner 
arborescence is a spanning tree that connects all points with M-paths to the origin. As in MMN, the 



aim is to find a minimum-weight network. For 2D-RSA, there is a polynomial-time approximation 
scheme (PTAS) [13] based on Arora's technique for approximating geometric optimization problems 
such as TSP pQ. It is not known whether 2D-MMN admits a PTAS. Arora's technique does not 
directly apply here as M-paths between terminals forbid detours and thus may not respect portals. 

1.3 Our contribution 

We first present a A(k — l)-approximation algorithm for the special case of 3D-MMN where the 
given terminals are contained in k > 2 planes parallel to the x—y plane; see Section [3] 

Our main result is an 0(n e )-approximation algorithm for cfD-MMN, for any e > 0. We first 
present the algorithm in detail for three dimensions; see Section [4} Since the algorithm for arbitrary 
dimensions is a straightforward generalization of the algorithm for 3D but less intuitive, we describe 
it in the appendix. 

Our 0(n £ )-approximation algorithm for ciD-MMN constitutes a significant improvement upon 
the best known ratio of 0(n 4 ' 5+£ ) for (general) directed Steiner forest [7j. We obtain this result 
by exploiting the geometric structure of the problem. To underline the relevance of our result, 
we remark that the bound of 0(n e ) is the best known result also for other directed Steiner- type 
problems such as DST [4| or even acyclic DST [18] . 

Our O (/^-approximation algorithm for the fc-planes case relies on recent work by Soto and 
Telha |17j . They show that, given a set of red and blue points in the plane, one can determine 
efficiently a minimum-cardinality set of points that together pierce all rectangles having a red point 
in the lower left corner and a blue point in the upper right corner. Combining this result with an 
approximation algorithm for 2D-MMN, yields an approximation algorithm for the 2-planes case. 
We show how to generalize this idea to k planes. 

2 Some Basic Observations 

We begin with some notation. Given a point p 6 I 3 , we denote the x-, y- and z-coordinate of p 
by x(p), y{p), and z(p), respectively. Given two points a and c in M 2 , let R(a,c) = {b € M 2 | 
x{a) < x(b) < x(c), y[a) < y(b) < y(c)} be the rectangle spanned by a and c. If a line segment is 
parallel to the x-, y-, or z-axis, we say that it is x-, y-, or z-aligned. In what follows, we consider 
the 3-dimensional case of the MMN problem, unless otherwise stated. 

2.1 Quadratic Lower Bound for Generating Sets in 3D 

Intuitively, what makes 3D-MMN more difficult than 2D-MMN is the following: in 2D, if the 
bounding box of terminals s and s' and the bounding box of t and t' cross (as in Fig. |lb[ ), then 
any s—s' M-path will intersect any t—tf M-path, which yields s-t' and t-s' M-paths for free (if s 
and t are the lower left corners of their respective boxes). A similar statement for 3D does not hold; 



M-paths can "miss" each other — even if their bounding cuboids cross; see Fig. Id 

Let us formalize this observation. Given a set T of terminals, a set Z of pairs of terminals is 
a generating set [11] if any network that M-connects the pairs in Z in fact M-connects all pairs 
of terminals. In 2D, any MMN instance has a generating set of linear size [TT]. Unfortunately 
this result does not extend to 3D. Below, we construct an instance that requires a generating set 
of size £2{n 2 ). The idea of using linear-size generating sets is exploited by several algorithms for 
2D-MMN |5|llj . The following theorem shows that these approaches do not easily carry over to 3D. 



Theorem 1. There exists an instance of 3D-MMN with n terminals that requires a generating set 
of size Q{n 2 ). 

Proof. We construct an instance that requires a generating set of size at least n 2 /4. The main idea 
of the construction is to ensure that n 2 /4 of the terminal pairs must use an edge segment unique 
to that specific pair. The input consists of two sets T and T", each with n/2 terminals, with the 
following coordinates: for < i < n/2, terminal ti G T is at (i, n/2 — i, n/2 — i) and terminal t[ G T" 
is at (n/2 + i,n — i,n — i). Figure ^1 shows the instance for n = 6. 




Fig. 2: The constructed network for n — 6. 

Consider any given generating set Z C T x T' such that there is a pair (t, t'), t G T and t' G T" 
that is not in Z. We now construct a specific network that contains M-paths between all terminal 
pairs in Z but no M-path between (t, t'). 

Consider any pair (t,t') G Z such that t = (i,j,k) G T and t' = (i',j',k r ) G T". The M-path 
from t to t' has three segments: an x-aligned segment from (i, j, k) to (i',j, k), a y-aligned segment 
from (i',j, k) to (i',f, k), and a z-aligned segment from (i',j', k) to (i',f, k'). To ensure an M-path 
between each generating pair ti,tj G T (similarly between t^t'- G X"), we add M-paths between 
each pair of consecutive terminals in T (similarly for T") as follows: we connect tj,tj+i G T by 
adding a ^-aligned segment from ti = (i,j, k) to (i,j, k — 1), a y-aligned segment to (i,j — l,k— 1), 
and an x-aligned segment to tj + i = (i + 1, j — 1, k — 1); see Fig. [2J 

It is easy to verify that, in this construction, the M-path between terminals t = (i,j,k) G T 
and t! = (i',j',k') G T" must use the y-aligned segment between (i',j,k) and (i',j',k). Since this 
segment is added only between terminal pairs that are present in the generating set Z, there is no 
M-path between terminals t G T and t' G T' which are not in Z. Thus, in order to obtain M-paths 
between all pairs of terminals in T U T', we need at least all of the n 2 /4 pairs in T x T'. □ 



2.2 Hanan Grid and Directional Subproblems 

First, we note that any instance of MMN has a solution that is contained in the Hanan grid, the 
grid induced by the terminals; see Fig. fl[a). Gudmundsson et al. [9] showed this for 2D; their proof 
generalizes to higher dimensions. In what follows, we restrict ourselves to finding feasible solutions 
that are contained in the Hanan grid. 



Second, to simplify our proofs, we consider the directional subproblem of 3D-MMN which 
consists of connecting all terminal pairs (£, t') such that t dominates t' , that is, x(t) < x(t'), 
y(t) < y(t'), z(t) < z(t'), and t / t'. We call such terminal pairs relevant. 

The idea behind our reduction to the directional subproblem is that any instance of 3D-MMN 
can be decomposed into four subproblems of this type. One may think of the above-defined di- 
rectional subproblem as connecting the terminals which are oriented in a north-east (NE) con- 
figuration in the x—y plane (with increasing z-coordinates) . Analogous subproblems exist for the 
directions NW, SE, and SW. Note that any terminal pair belongs to one of these four categories 
(if seen from the terminal with smaller z-coordinate). 

The decomposition extends to higher dimensions d, by fixing the relationship between (t, t') for 
one dimension (for example, z), and enumerating over all possible relationships for the remaining 
d — 1 dimensions. This decomposes <iD-MMN into 2 subproblems, which is a constant number 
of subproblems as we consider d to be a fixed constant. 

This means that any /^-approximation algorithm for the directional subproblem leads to an 
0(p)-approximation algorithm for the general case. Thus we can focus on designing algorithms for 
the directional subproblem. 

Observation 1 Any instance of MMN can be decomposed into a constant number of directional 
subproblems. Thus a p- approximation algorithm for the directional subproblem leads to an O(p)- 
approximation algorithm for MMN. 

2.3 Relation to Steiner Problems 

We next show that there is an approximation-preserving reduction from directional 3D-MMN to 
the directed Steiner forest (DSF) problem, which by Observation [II carries over up to a constant 
factor, to general 3D-MMN. 

Let T be a set of n points in M 3 . Let H be the Hanan grid induced by T. We consider H as an 
undirected graph where the length of each edge equals the Euclidean distance between its endpoints. 
We orient each edge in H so that, for any edge (p,p') in the resulting digraph H', the start node p 
dominates the end node p' . We call H 1 the oriented Hanan grid of T. Now let (£, t') be a relevant 
pair of points in T, that is, t dominates t' . Any M-path in H connecting t to t' corresponds to a 
directed path in H' from t to t' . The converse also holds: every directed path in H' corresponds to 
an M-path in H. 

Let / be an instance of directional 3D-MMN and let I' be an instance of DSF where the input 
graph is H' and where every relevant terminal pair of / has to be connected. Then, each feasible 
solution N of I contained in H corresponds to a sub-graph N' of H' that connects every relevant 
terminal pair, and is therefore a feasible solution to I'. It is easy to see that N' has the same cost as 
N, as N' uses the oriented version of each edge of N. Conversely, every feasible solution N' for I' 
corresponds to a subgraph N of H that M-connects every relevant terminal pair. Therefore, A is a 
feasible solution to I' with the same cost as A'. This establishes an efficiently computable one-to- 
one correspondence between feasible solutions to I that are contained in H and feasible solutions to 
I'. Since there is an optimum solution to I contained in H [9j, this is an approximation-preserving 
reduction from directional 3D-MMN to DSF. 

By means of the above transformation of the Hanan grid into a digraph, we also obtain an 
approximation-preserving reduction from 3D-RSA to DST. We use this later in Section [4] to develop 
an approximation algorithm for 3D-MMN. Let I be an instance of 3D-RSA given by a set T of 



terminals with non-negative coordinates that are to be M-connected to the origin o. We construct 
an instance I' of DSF as above where {o} x T is the set of node pairs to be connected. Note that 
any feasible solution to I' is, without loss of generality, a tree. Hence, I' is an instance of DST with 
root o. All in all, we have an approximation-preserving reduction from 3D-RSA to DST. 

3 The fc-Plane Case 

In this section we consider 3D-MMN, under the assumption that the set T of terminals is contained 
in the union of k > 2 planes E\, . . . , Ek that are parallel to the x-y plane. Of course, this assumption 
always holds for some k < n. We present a 4(/c — l)-approximation algorithm, which outperforms 
our algorithm for the general case in Section^ if k £ o(n 6 ). 

Let iV pt be some fixed minimum M-network for T, let N^ be the set of all z-aligned and all 
y-aligned segments in -/Vop t , and let N^S be the set of all z-aligned segments in -/V pt. Let OPT 
denote the weight of -/Vo P t • Clearly, OPT does not depend on the specific choice of iV pt ; the weights 
of Nq°1 and N™1, however, may depend on iV op t. For i € {1, . . . , k}, let Tj = T fl Ei be the set of 
terminals in plane Ei. Further, let T xy be the projection of T onto the x~y plane. 

Our algorithm consists of two phases. Phase I computes a set N of horizontal (that is, x- and 
y-aligned) line segments, phase II computes a set iV ver of vertical (that is, z-aligned) line segments. 
Finally, the algorithm returns the set N = N hm U N ver . 

Phase I is simple; we compute a 2-approximate M-network N xy for T xy (using the algorithm of 
Guo et al. [10]) and project N xy onto each of the planes E\, . . . , E^. Let N hor be the union of these 
projections. Note that JV M-connects any pair of terminals that lie in the same plane. 

Observation 2 \\N hm \\ < 2/c||7V h ° t r ||. 

Proof. The projection of N^°l to the x-y plane is an M-network for T xy . Hence, H-ZV^H < 2A^° t r . 
Adding up over the k planes yields the claim. □ 

In Phase II, we construct a pillar network by computing a set iV ver of vertical line segments, so-called 
pillars, of total cost at most 4(fc — 1 ) 1 1 A^T J; 1 1 . This yields an overall approximation factor of 4(k — 1) 
since ||iV hOT U N™\\ < 2k\\N^\\ + 4(fc - l)||iV V p C t r || < 4(fc - l)(||iV h p-|| + ||iV v p e t l) < A(k - l)OPT. 

Below we describe Phase II of our algorithm for the directional subproblem that runs in direction 
north-east (NE) in the x-y plane (with increasing z-coordinates) . For this directional subproblem, 
we construct a pillar network iVJ? r of weight at most (k — l)||^V"optll that, together with ./V" hor , M- 
connects all relevant pairs. We solve the analogous subproblems for the directions NW, SE, and SW 
in the same fashion. Then JV ver is the union of the four partial solutions and has weight at most 
4(k-l)\\N™l\\, as desired. 

Our directional subproblem is closely linked to the (directional) bichromatic rectangle piercing 
problem (BRP), which is defined as follows. Let R and B be sets of red and blue points in M 2 , 
respectively, and let 1Z(R, B) denote the set of axis-aligned rectangles each of which is spanned by 
a red point in its SW-corner and a blue point in its NE-corner. Then the aim of BRP is to find a 
minimum-cardinality set PcR 2 such that every rectangle in 1Z(R, B) is pierced, that is, contains 
at least one point in P. The points in P are called piercing points. 

The problem dual to BRP is the (directional) bichromatic independent set of rectangles problem 
(BIS) where the goal is to find the maximum number of pairwise disjoint rectangles in 1Z{R,B), 
given the sets R and B. 
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Recently, Soto and Telha [17] proved a beautiful min-max theorem saying that, for 1Z(R, B), the 
minimum number of piercing points always equals the maximum number of independent rectangles. 
This enabled them to give efficient exact algorithms for BRP and BIS running in 0(n 2,5 ) worst-case 
time or Ofa"') expected time, where the O-notation ignores polylogarithmic factors, 7 < 2.4 is the 
exponent for fast matrix multiplication, and n = \R\ + \B\ is the input size. 



The details of Phase II appear, for k = 2 planes, in Section 3.1 and, for k > 2 planes, in 
Section 3.2 Algorithm [T] summarizes of our /c-planes algorithm. 



Algorithm 1 /c-Planes Algorithm 



Input: Set T of terminals contained in the union of planes E\, . . . ,Ek, all parallel to the x-y plane. 

1: Let T xy be the projection of T onto the x-y plane 
2: Phase I: 

Compute N xy , a 2-approximate M-network for T xy using the algorithm of Guo et al. [10] . 

Let N hoT be the union of the projections of N xy onto each of the planes E\ , . . . , E^ . 
3: Phase II: 

If k = 2, construct a pillar network JV vcr by Algorithm [2] see Section [3. 1| 

Otherwise, construct a pillar network iV vor by Algorithm [3] see Section [3.2| 
Output: iV hor UAr vcr . 



3.1 Pillar Network for Two Planes 

Our phase-II algorithm for two planes is very simple. We sketch it first in order to provide some 
intuition for the fc-planes case. Let the terminals in T± be red and those in T2 be blue. Ignore the 
^-coordinates of the terminals. Then the relevant red-blue point pairs span exactly the rectangles 
in TZ{T\,T2), which we call relevant, too. 

Algorithm 2 Pillar network of the directional subproblem for k = 2 planes 
Input: Sets Ti C -Bi and T2 C -E2 of terminals. 

1: Color Ti red and T 2 blue. 

2: Ignoring z-coordinates of terminals, let TZ(Ti,T2) be the set of rectangles spanned by relevant red-blue pairs. 

3: Compute a minimum piercing P of TZ(Ti , T2) such that for each relevant red-blue pair (r, b) G Ti x T2 the piercing 
point for (r, 6) lies on an r-b M-path in N xy , as described in Lemma [5] 

4: Erect pillars from E\ to E2 at each piercing point p £ P; let N^ be the resulting set of pillars. 
Output: AX. r . 

Our algorithm (Algorithmic]) consists of two steps. First, we compute a minimum piercing P of 
1Z{Ti,T2) using the algorithm of Soto and Telha [17] . Second, we move each piercing point p G P 
to a new position p — a nearby junction of N xy — and erect, at p, a pillar connecting the two planes. 
Let P be the set of piercing points after the move, and let N^ be the corresponding set of pillars. 

Lemma 1. It holds that \\N$g\\ < \\N^\\. 

Proof. It is easy to see that |P| = \P\. Integrating over the distance d of the two planes yields 
||AJP r || = \P\ • d = \P\ • d < ||A~op{:||. The last inequality is due to the fact that P is a minimum 
piercing of 1Z(T\,T2) and that the pillars in N™l pierce TZ(Ti,T2) — otherwise A" op t would not be 
feasible. □ 



Now we turn to feasibility. We first detail how we move each piercing point p to its new position p. 
For the sake of brevity, we identify terminals with their projections to the x—y plane. Our description 
assumes that we have at our disposal some network M (such as N xy ) connecting the relevant pairs 
in T xy . 

For a piercing point p E P, let A p be the intersection of the relevant rectangles pierced by p; see 
Fig. [31 Clearly, p e A p . Note that the bottom and left sides of A p are determined by terminals £\y 
and is to the west and south of A p , respectively. Symmetrically, the top and right sides of A p are 
determined by terminals £e and £n to the east and north of A p , respectively. Terminals £\y and ts 
may coincide, and so may £e and £n- It is easy to see that the network M contains an M-path 7Tsn 
connecting t$ and % and an M-path 7Twe connecting tw and £e- The path 7Tsn goes through the 
bottom and top sides of A p and 7Twe goes through the left and right sides. Hence, the two paths 
intersect in a point p £ A p . This is where we move the original piercing point p. 




Fig. 3: Paths 7tsn and 7twe meet in a point p in A p . 



Since p E A p , the point p pierces the same relevant rectangles as p, and the set P = {p \ p £ P} 
is a (minimum) piercing for the set of relevant rectangles. 

Lemma 2. Let 1Z(R, B) be an instance of BRP and let M be a network that M-connects every 
relevant red-blue point pair. Then we can efficiently compute a minimum piercing of 1Z(R, B) such 
that M contains, for every relevant red-blue point pair (r, b) in Rx B, anr-b M-path that contains 
a piercing point. 

Proof. We use the algorithm of Soto and Telha [TT] to compute a minimum piercing P of 1Z(R, B). 
Then, as we have seen above, P is a minimum piercing of 1Z(R, B), too. Now let (r, b) be a relevant 
red-blue pair in R x B, and let p £ P be a point that pierces R(r, b). Clearly, p pierces R(r, b), too. 
As we have observed before, both p and p lie in A p . 

Since (r, b) is a relevant pair, r lies to the SW of A p and b to the NE; see Fig. 4a We prove 
that M contains an r—p M-path; a symmetric argument proves that M also contains a p-b M-path. 
Concatenating these two M-paths yields the desired r-b M-path since r lies to the SW of p and p 
lies to the SW of b. Recall that p lies on the intersection of the tw~^E M-path 7T\ve and the £s - ^N 
M-path 7TgN, where t-^f, £e, *s> *n are the terminals that determine the extensions of A p ; see Fig.|3j 
To show that M M-connects r and p, we consider two cases. 



Case I: r E R(tw, ts); see Fig. 4a According to our assumption, M contains some r-b M-path n. 
Then tt must intersect 7Twe or 7Tsn at some point x to the SW of p. Thus, we can go, in a monotone 
fashion, along tt from r to x and then along 7Twe or 7Tsn from x to p. This is the desired r—p M-path. 



Case II: r lies to the SW of £w ° r fej see Fig- 4b In this case M contains M-paths from r to £\y 
and to is- If r lies to the SW of tyy, we can go, again in a monotone fashion, from r to fr\y and 





°b 



(a) Case I: we can go from r via x 
to p. 



(b) Case II: we can go from r via iw or is 
to p. 



Fig. 4: Sketches for the proof of Lemma 2 

then along 7T\ve from £w to p. Otherwise, if r lies to the SW of ts, we can go from r to ts an d then 
on 7Tsn from is to p. 

Since these are the only two possibilities, this concludes the proof. □ 

Lemmas [I] and ^1 (with R = T\, B = T2, and M = N xy ) yield the following. 
Theorem 2. We can efficiently compute a 4- approximation for the 2-plane case. 



3.2 Pillar Network for k Planes 

Now we show how our phase-II algorithm generalizes to k planes. As in the 2-planes case, we 
restrict ourselves to the directional subproblem and construct a pillar network NJf* of weight 
at most (k — l)\\N™ T t \\. As we have argued at the beginning of Section [3l this suffices to prove 
Theorem [3l 

Theorem 3. There exists a 4(fc — 1)- approximation algorithm for 3D-MMN where the terminals 
lie in the union of k >2 planes parallel to the x-y plane. 

Our pillar-placement algorithm (Algorithm k}| is as follows. Let i G {1, . . . , k — 1}. We construct 
an instance 2j of BRP where we two-color T xy such that each point corresponding to a terminal 
of some plane Ej with j < i is colored red and each point corresponding to a terminal of some 
plane Ej/ with j' > i + 1 is colored blue. For Xj, we compute a minimum piercing Pi according to 
Lemma I2I with M = N xy . In other words, for any relevant pair (tj,tji), there is some M-path in 
N xy that contains a piercing point of Pi. We choose i k G {1, . . . , k — 1} such that Pi* has minimum 
cardinality. This is crucial for our analysis. At the piercing points of P^, we erect pillars spanning 
all planes E\, . . . , E^. Let Ni* be the set of these pillars. We now show that iVj*, along with N , 
creates a feasible network for any relevant terminal pair (tj,t'A such that j < i* and j' > i* + 1. 

Lemma 3. The network N hov U Ni* M-connects any relevant terminal pair in Tj x T,y with j < i* 
and j' > i* + 1. 

Proof. Consider a pair (tj,tj') in Tj x Tji as in the statement. We construct an M-path from tj 
to tj' as follows. We know that there exists an M-path ir that connects the projections of tj and ty 
in N xy and contains a piercing point p of Pi* . Therefore, we can start at tj and follow the projection 
of 7r onto plane Ej until we arrive at p. Then we use the corresponding pillar in Ni* to reach the 
plane Eji, where we follow the projection of 7r (onto that plane) until we reach tj'. □ 



Algorithm 3 Pillar network for the directional subproblem for k > 2 planes 

Input: Sets T s C E s , . . . , T t C E t of terminals with s < t (initially s = 1 and t — k). 

1: Let T' be the projection of T s U ■ ■ • U T t onto the x-y plane. 

2: for each i G {s, . . . , i} do 

3: Let Xi be an instance of BRP where each point in T" , corresponding to a terminal in Tj with j < i, is colored 
red and each point in T' , corresponding to a terminal in Tji with j' > i + 1, is colored blue. 

4: Compute a minimum piercing Pi according to Lemma |2| with M — N xy . 

5: end for 

6: Choose i* £ {s, . . . ,£} such that Pi* has minimum cardinality. 

7: Let Ni* be the set of pillars erected at each piercing point of Pi* , spanning planes E 3 , . . . ,Et. 

8: Let N<i* be the output of this algorithm applied recursively to T s , . . . , Ti* . 

9: Let N>i* be the output of this algorithm applied recursively to Ti*+i, . . . ,Tt. 
Output: Ni* U N<i* U N>i* 



In order to also M-connect relevant terminal pairs in Tj x Tji, where either (j < i* and j' < i*) 
or (j > i*-\-l and j' > i* + l), we simply apply the pillar-placement algorithm recursively to the sets 
Ti, ... ,Tj* and Tj* + i, . . . ,T^. This yields the desired pillar network iVJ? r . By Lemmapl AQ^f U iV hor 
is feasible. Next, we bound ||iVj*||. 

Lemma 4. Let M be an arbitrary directional Manhattan network for T , and let M vev be the set of 
vertical segments in M. Then the pillar network Ni* has weight at most ||M ver ||. 

Proof. Without loss of generality, we assume that M is a subnetwork of the Hanan grid [9j . We may 
also assume that any segment of M ver spans only consecutive planes. For 1 < i < j < k, let Mij 
denote the subnetwork of M vcr lying between planes Ei and Ej. Let ckj be the vertical distance 
between planes Ei and Ej. 

We start with the observation that, for any j = 1, . . . , k — 1, the network Mjj + \ is a set of 
pillars that forms a valid piercing of the piercing instance Xj (defined right after Theorem ^). 
Hence, iMjj+i) > \Pj\ > |Pj*|, which implies the claim of the lemma as follows: 

fe-i fc-i fc-i 

||M ver || = J^ ||M iii+1 || = Y, \M jJ+ i\ ■ d h3+1 > Y, |iV| • d SJ+ i = \P t *\ ■ di, fc = ||P*||. 

D 

It is crucial for our construction that the pillars constructed recursively span either E\ , . . . , Ei* 
or Ei*+i, . . . ,Ek, but not all planes. For 1 < j < j' < k, let weight z (j, j') denote the weight of 
the vertical part of the network produced by the above pillar-placement algorithm, when applied 
to planes Ej, . . . ,Eji recursively. For technical reasons we set weight z (j,j) = 0. Now assume that 
j < j' and that the algorithm makes the partition at plane Ey with j < i 1 < j' when planes 
Ej, . . . , Ej/ are processed. By means of Lemma |4l we derive the recursion 

weighty, f) < \\Mjj, || + weighty, i') + weight^' + 1, j') , (1) 

which holds for any M-network M for T. We now claim that 

weight z (j,j')<(j'-j)\\M jd ,\\. 

10 



Our proof is by induction on the number of planes processed by the algorithm. By the induc- 
tive hypothesis, we have that weighty (j,i') < (i' — j)\\Mj t ii\\ and weight 2 (i' + l,f) < (j' — i' — 
l)||Mj/ +1 j>||. We plug these expressions into the recursion [l] Since ||M^/|| + ||Mj/_|_ lj /|| < ||Mjj/|| 
and weight z (Z, I) = for any I £ {1, . . . , k}, the claim follows. 

We conclude that the weight of the solution produced by the algorithm, when applied to all 
planes E u ...,E k ,is bounded by weighty (1, k) < (k - l)||Mi )fe [| = (k - l)||M vor ||. This completes 
the proof of Theorem [3| 

4 The General Case 

In this section, we present an approximation algorithm, which we call the grid algorithm, for the 
general 3D-MMN problem. Our main result is the following. 

Theorem 4. For any e > 0, there exists an 0{n £ )- approximation algorithm for 3D-MMN that, 
given a set of n terminals, runs in n°(V e ) time. 

This result is better than the one in the previous section if the given set of terminals is distributed 
over uo(n £ ) horizontal planes. Moreover, the approach in this section extends to higher dimensions; 
see appendix. 

For technical reasons, we assume that the terminals are in general position, that is, any two 
terminals differ in all three coordinates. By Observation [T] it suffices to describe and analyze the 
algorithm for the directional subproblem. 

4.1 The 3D Grid Algorithm 

We begin the description with a high-level summary. To solve the directional subproblem, we place 



a 3D grid that partitions the instance into a constant number of cuboids; see Fig. 5a Cuboids that 
differ in only two coordinates form slabs. We connect terminals from different slabs by M-connecting 
each terminal to the corners of its cuboid and by using the edges of the grid to connect the corners. 
We connect terminals from the same slab by recursively applying our algorithm to the slabs. 

Step 1: Partitioning into cuboids and slabs. Consider the bounding cuboid C of T and set c = 3 l < £ . 
Partition C by 3(c — 1) separating planes into c x c x c axis-aligned subcuboids C^k with i, j, k E 
{1, . . . ,c}. The indices are such that larger indices mean larger coordinates. Place the separating 
planes such that the number of terminals between two consecutive planes is at most n/c. This can be 
accomplished by executing a simple plane-sweep for each direction x, y, z, and by placing separating 
planes after every n/c terminals. Here we exploit our general-position assumption. The edges of 
the resulting subcuboids — except the edges on the boundary of C, which we do not need — induce 
a three-dimensional grid Q of axis-aligned line segments. We insert Q into the solution. 

For each i G {1, . . . , c}, define the x-aligned slab, Cf , to be the union of all cuboids Cy& with 
j, k £ {1, . . . , c}. Define y-aligned and z-aligned slabs C y -, Cf. analogously; see Fig. 



5b 



Step 2: Add M-paths between different slabs. Consider two cuboids C^k and Cityy with i < i! , 
j < j', and k < k! . Any terminal pair (t, t') £ C^ X Ci'jiy can be M-connected using the edges 



of Q as long as t and t' are connected to the appropriate corners of their cuboids; see Fig. 5c To 
this end, we use the following patching procedure. 

Call a cuboid C^k relevant if there is a non-empty cuboid Ci'j'k' with i < i! , j < f , and k < k'. 
For each relevant cuboid Cijk, let pijk denote a corner that is dominated by all terminals inside C^k- 
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(a) Example of a grid with two (b) Examples of X-, y-, and z-aligned ( c ) M-paths using patching 

cuboids from different slabs. slabs. anc l n-j-jj edges. 

Fig. 5: Illustrations for the grid algorithm: cuboids, slabs, and patches. 



We define up-patching Cijk to mean M-connecting every terminal in Cijk to pijk ■ We up-patch Cijk 
by solving (approximately) an instance of 3D-RSA with the terminals in C^k as points and pijk as 
origin. We define down-patching analogously; cuboid C^k is relevant if there is a non-empty cuboid 
Ci'j'k' with i > i! , j > j' , k > k'; we let pijk be the corner that dominates all terminals in Cijk- 

We complete this step by inserting the up-patches and the down-patches of all relevant cuboids 
into the solution. 

Step 3: Add M-paths within slabs. To M-connect relevant terminal pairs that lie in the same slab, we 
apply the grid algorithm (steps 1-3) recursively to each slab Cf, Cj, and C%. with i,j, k G {1, . . . , c}. 

4.2 Analysis 



We first show that the output of the algorithm presented in Section 4.1 is feasible, then we establish 
its approximation ratio of 0(n e ) and its running time of n ^ 1 ' 6 ' for any e > 0. In this section, OPT 
denotes the weight of a minimum M- network (not the cost of an optimal solution to the directional 
subproblem) . 

Lemma 5 (Feasibility). The grid algorithm M-connects all relevant terminal pairs. 

Proof. Let (t, t') be a relevant terminal pair. First, suppose that t and t' lie in cuboids of different 
slabs. Thus, there are i < i',j < j',k < k! such that t G Cijk and t' £ Ciijiy. Furthermore, Cijk 
and Ciij/y are relevant for up- and down-patching, respectively. When up-patching, we solve an 
instance of RSA connecting all terminals in C^k to pijk- Similarly, down-patching M-connects t' 
to pi'j'k'- The claim follows as Q M-connects pijk and pi'j'k'- 

Now, suppose that t and t' lie in the same slab. As the algorithm is applied recursively to each 
slab, there will be a recursion step where t and t' lie in cuboids in different slabs. Here, we need 
our general-position assumption. Applying the argument above to that particular recursive step 
completes the proof. □ 

Approximation ratio. Next, we turn to the performance of our algorithm. Let r(n) be its approxi- 
mation ratio, where n is the number of terminals in T. The total weight of the output is the sum 
of ||£/||, the cost of patching, and the cost for the recursive treatment of the slabs. We analyze each 
of the three costs separately. 
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The grid Q consists of all edges induced by the c 3 subcuboids except the edges on the boundary 
of C. Let £ denote the length of the longest side of C. The weight of Q is at most 3(c — l) 2 ^, which 
is bounded by 3c 2 OPT as I < OPT. 

Let r pa tch(^) denote the cost of patching all relevant cuboids in step 2. Lemma M (given below) 
proves that r pa t c h(^) = 0(n 6 ) OPT. 

Now consider the recursive application of the algorithm to all slabs. Recall that A op t is a fixed 
minimum M- network for T. For i £ 1, . . . , c, let OPTf be the optimum cost for M-connecting all 
(not only relevant) terminal pairs in slab Cf. Define OPT^ and OPTf analogously. 

Slightly abusing of notation, we write A op t n Cf for the set {sCiCf \ s £ N opt } of line segments 
of iV pt intersected with slab Cf. Observe that N op t n Cf forms a feasible solution for Cf. Thus, 
OPTf < H-ZVopt n Cf\\. By construction, any slab contains at most n/c terminals. Hence, the total 
cost of the solutions for slabs Cf , . . . , C% is at most 



C C 



J2r{- c )OPTf<r^)j:\\N opt nCn<r^)OFT 

i=l i=l 

Clearly, the solutions for the y- and z-slabs have the same bound. 

Summing up all three types of costs, we obtain the recursive equation 



r(n)OPT < 3cWT + r patch (ra)OPT + 3r I -J OPT . 

Hence, r(n) = 0(n max ^ e ' gc3 J). Plugging in c = 3 1 ' 6 yields r(n) = 0(n e ), which proves the approx- 
imation ratio claimed in Theorem |U 

Lemma 6. Patching all relevant cuboids costs r pa tch(^) £ 0(n £ ) OPT. 

Proof. First note that it suffices to consider up-patching; the down-patching case can be argued 
analogously. 

Lemma [7] shows the existence of a near-optimal M-network that up-patches all relevant cuboids. 
Lemma [8] shows that by reducing the patching problem to 3D-RSA, we can find such a network of 
cost 0(p)OPT, where p is the approximation factor of 3D-RSA. 



We argued in Section 2.3 that there exists an approximation-preserving reduction from 3D-RSA 
to DST. DST, in turn, admits an 0{n e )-approximation for any e > j3]- Hence, the cost of up- 
patching is indeed bounded by 0(n £ )OPT. □ 

We now turn to the two lemmas that we just used in the proof of Lemma [6| For our analysis, 
we need the network N' that is the union of Q with A^pt and the projections of iV pt onto every 
separating plane of Q. Since there are 3(c — 1) separating planes and, as we have seen above, 
||a|| < 3c 2 OPT, it holds that \\N'\\ < 3(c 2 + c)OPT = O(OPT). 

Lemma 7. There exists an M-network of total cost at most 3(c 2 + c)OPT that up-patches all 
relevant cuboids. 

Proof. We claim that N' up-patches all relevant cuboids. To this end, let t £ Cijk and let t' G Cy^y 
with i < i',j < f, k < k' . Follow the M-path connecting t and t', starting from t. This path must 
leave C^k at a certain point p, which lies on some face F of C^k- Face F, in turn, lies on some 
separating plane S of the grid Q. From now on follow the projection of the M-path from p to t' on 
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plane S. This projected path must leave the face F, since t' lies in Ci'j>k> with i < i' ,j < j', k < k', 
and the projection of t' onto S must therefore lie outside of F. Moreover, the point p' where this 
path leaves F must lie on an edge of Cijk incident to pijk- Hence, we obtain a t-pijk M-path by 
going from t to p, from p to p' and then from p' to Pijk- 

Lemma 8. Given a number p > 1 and an efficient p- approximation of 3D-RSA, we can efficiently 
up-patch all relevant cuboids at cost no more than 12(c 2 + c)pOPT. 

Proof. In Lemma [7| we showed the existence of a network N' that up-patches all relevant cuboids 
at low cost. Now consider an arbitrary relevant cuboid C^k- Clearly N' n C^k up-patches C^k- 
Hence OPT" 1 ^, < \\N' (1 Cijk\\, where OPTS, denotes the cost of a minimum up-patching of Cijk- 
The problem of optimally up-patching C^k is just an instance Jyfe of 3D-RSA in which all terminals 
in C^k have to be connected by an M-path to pijk- Applying the factor-/) approximation algorithm 
for 3D-RSA to each instance Iijk with C^k relevant, we patch at total cost at most 

p J2 0PT ^ ^ P E W N 'nC iik \\ < 4p||iV'||. 

Cijk relevant Cijk relevant 

The last inequality follows from the fact that each edge of N' occurs in at most four cuboids. The 
lemma follows since [|JV'|| < 3(c 2 + c)OPT. D 

Running time. Finally, we analyze the running time. Let T{n) denote the running time of the 
algorithm applied to a set of n terminals. The running time is dominated by patching and the 
recursive slab treatment. Using the DST algorithm of Charikar et al. [1], patching cuboid C% requires 
time n i , where n, is the number of terminals in Cj. As each cuboid is patched at most twice 
and there are c 3 cuboids, patching takes 0(c 3 )n ^' £ ' = n \ l e ' time. The algorithm is applied 
recursively to 3c slabs. This yields the recurrence T{n) = 3cT(n/c) + n °' 1 / £ J ) which leads to the 
claimed running time. 

This completes the proof of Theorem |4j 

5 Open Problems 

We have presented, for any e > 0, a grid-based 0(n £ )-approximation algorithm for dD-MMN. 
This is a significant improvement over the ratio of 0(n 4 ' 5+e ) which is achieved by reducing the 
problem to DSF. For 3D, we have described a 4(A; — l)-approximation algorithm for the case when 
the terminals lie on k > 2 horizontal planes. This outperforms our grid-based algorithm when 
k G o{n £ ). Whereas 2D-MMN admits a 2-approximation |5)10|15j . it remains open whether 0(1)- 
or 0(logn)-approximation algorithms exist for higher dimensions. 

Our 0(n £ )-approximation algorithm for dD-MMN solves instances of (iD-RSA for the subprob- 
lem of patching. We conjecture that (iD-RSA admits better approximation ratios. While this is an 
interesting open question, a positive result would still not be enough to improve our approximation 
ratio, which is dominated by the cost of finding M-paths inside slabs. 

The complexity of the undirectional bichromatic rectangle piercing problem (see Section pi) is 
still unknown. Currently, the best approximation has a ratio of 4, which is (trivially) implied by the 
result of Soto and Telha [17J. Any progress would immediately improve the approximation ratio of 
our algorithm for the &-plane case of 3D-MMN (for any k > 2). 
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Appendix: Extension to Higher Dimensions 

We now describe the approximation algorithm for <iD-MMN, for d > 3, as a generalization of the 
3D-MMN idea from Section H 

Theorem 5. For any fixed dimension d and for any e > 0, there exists an 0(n £ ) -approximation 
algorithm for dD-MMN. 

Large parts of the algorithm and the analysis are straightforward generalizations of the algorithm 



for 3D-MMN. The presentation of both follows closely the 3D case. However, Lemma 11 where the 
cost of the patching procedure is bounded, requires non-trivial additional insights. 

As in the 3D case we decompose the overall problem into a constant number of instances of 
the directional subproblem. The directional subproblem consists in M-connecting any terminal pair 
(t, t') such that x l {t) < x l (t') for any i G {1, . . . , d}. Here, we use x l (p) to denote the i-th coordinate 
of point p G M. d . We can decompose the general problem into 2 d_1 directional subproblems. Once 
again, we assume that the terminals are in general position. 

The dD Grid Algorithm 

We begin the description with a high-level summary. To solve the directional dD-MMN problem 
we place a dD-grid which partitions the instance into cuboids and slabs. Terminal pairs lying in 
different slabs are handled by M-connecting each terminal to the corner of its cuboid and then using 
the edges of grid. Terminal pairs from the same slab are M-connected by applying the algorithm 
recursively to all slabs. Each slab contains only a constant fraction of the terminals. 
Step 1: Partitioning into cuboids and slabs. Consider the bounding cuboid C for the set T of 

terminals and choose a large constant c = d > E . For each dimension i G {1, . . . , d} we choose c + 1 
separating planes determined by values x\ < ■■• < x l c+1 . Planes x\ and x l c+1 coincide with the 
boundary of C in dimension i. The separating planes for dimension i partition C into c slabs C*-, 

where j G {1, . . . , c}. Slab C* is the set of points p G C such that x\ < Xi(p) < x\ . We place the 
separating planes so that each slab contains at most n/c terminals. Altogether we have d(c + 1) 
separating planes. 

Let ji, . . . ,jd G {1, ... , c}. The subcuboids Cj lj ,,,j d is the set of points p G C such that x l - < 
x l (p) < x l j_ +l for each i G {1, . . . , d}. 

Consider the i th dimension, i G {1, . . . , d} and integers jk G {2, . . . , c} for each k G {1, . . . , d} — 
{i}. Let s be the axis-parallel line segment that contains all points p£C such that x k (p) = Xj k for 
each k ^ i. We call s a grid segment for dimension i. The grid Q is the set of all grid segments and 
there are d(c — l) d ~ 1 grid segments in Q. 

Given jk G {2, . . . , c} for each k G {1, . . . , d} we call (xj a , . . . , x d , d ) a grid point of Q and there 
are (c — 1) grid points in total. 

Step 2: Add M-paths between different slabs. Consider two cuboids Cj lt ,,,j d and Cj',...,j' with ji < j\ 
for each i G {1, . . . ,d}. Any pair of terminals t G Cy ,-' and t! G C,-/ ,-' can be M-connected 
using the segments of Q as long as t and t are suitably connected to the corners (grid points) of 
their cuboids. We use patching (described below) to connect all terminal to the corners of their 
cuboid. 

Patching: Call a cuboid Cj ly ...j d relevant if there is a cuboid Cji y that contains at least one 
terminal and satisfies ji < j[ for each % G {1, . . . , d}. For each relevant cuboid Cj lt ,,,j d , let Pj 1> ...,j d 
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denote the grid point (x*- +1 , . . . ,xq +1 ). Up-patching Cj lt ,,,j d means to M-connect every terminal 
in i € Cji,...,id to Pj 1 ,...,j d - Up-patch Cj lt .„j d by solving the gED-RSA problem with the terminals 
inside Cj t j d as the terminals and Pj 1 ,... 7 j d as the origin. 

Down-patching is defined analogously; cuboid Cj lt ...j d is relevant if there is a non-empty cuboid 
Cj> t ji with j[ < ji, i = 1, . . . , d and using grid point p', • := (x* i+1 , . . . , x^ +1 ) as origin instead 

The output of this step is the union of grid Q with a network that up-patches and down-patches 
all relevant cuboids. This produces M-paths between all terminal pairs in different slabs. 

Step 3: Add M-paths within slabs. To also connect terminal pairs that lie in a common slab we apply 
the algorithm (Steps 1-3) recursively to each slab Cj with i G {1, . . . , d} and any j G {1, . . . , c}. 

Analysis 

We now show that the algorithm presented above yields a feasible solution to directional dD-MMN, 
with cost at most 0(n e )OPT, for any e > 0. Here, OPT denotes the cost of an optimum solution to 
the general dD-MMN instance rather than the minimum cost OPT' achievable for the directional 
subproblem. The reason is that the cost of the grid Q is generally not related to OPT' but to OPT. 
We finish the section by arguing that the running time of the algorithm is n°( l ' £ > . 

Lemma 9 (Feasibility). The dD grid algorithm M-connects all relevant terminal pairs. 

Proof. Let (t,t') be a relevant terminal pair. First suppose that t £ Cj lt ,„j d and t' S Cy ,■/ where 
ji < j[ for all i 6 {1, . . . ,d}. Hence, Cj lt ...j d and Cj' j' are relevant for up-patching and down- 
patching, respectively. Consider the corners pj, ,• . of C,-, , , and the corner p'., ., of C,-/ ,-' . 
In the up-patching step of our algorithm we solve an RSA problem with terminals of Cj ll ..,j d as 
the input points, and corner Pj 1 ,....j d as the origin. By definition, an RSA solution M-connects t to 
Pji,—,jd- Similarly, down-patching M-connects t' to p'-, ., . It follows that t and t! are connected, 
since Pj lt ...,j d and p'., ., are M-connected via grid Q. Both terminals are even M-connected since 
additionally t < Pj 1 ....,j d < p'i ■> < t', where < denotes the domination relation between points. 

Now suppose t and t' lie in the same slab. As the algorithm is applied recursively to each slab 
there will be a recursion step where t and t' will lie in cuboids in different slabs. Here, we need 
our assumption of general position. Applying the argument above to that particular recursive step 
completes the proof. □ 

Approximation ratio. Let r{n) denote the approximation ratio of our algorithm where n is the 
number of terminals in T. The total cost of our solution consists of the cost for the grid Q, the cost 
of up-patching and down-patching all relevant cuboids, and the cost for the recursive treatment of 
the slabs in all d dimensions. We analyze each of these costs separately. 

The grid Q consists of the d(c— l) d ~ l grid segments. The length of any grid segment s is a lower 
bound on OPT. This holds because there are two terminals on the boundary of C whose Li-distance 
is at least the length of s. It follows that the cost of the grid is bounded by d(c — l) d_1 OPT. 



Let r pa tch(ra) denote the cost of patching all relevant cuboids as is done in Step 2. Lemma 10 
(given below) proves that r patc h(n) = 0(n £ )OPT. 

Now consider the recursive application of the algorithm to all slabs C!j, where i £ {1, . . . ,d} 
and j £ {1, . . . ,c}. First recall that we placed the separating planes so that \C l A < n/c for any 
i 6 {1, . . . , d} and any j £ {1, . . . , c}. 
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Consider dimension i G {1, . . . ,d}. Let OPT*- be the optimum cost for M-connecting all (not 
only relevant) terminal pairs in slab CL where j G {1, . . . , c}. Slightly abusing notation, we write 
-Wopt H Cj for the set of line segments of iV pt that are completely contained in the slab C* . Observe 
that iVopt H Cj, forms a feasible solution for Cj. Thus OPT*- < ||JV pt n Cj||. Each such Cj contains 
at most n/c terminals, and therefore the total cost of the solutions for the all slabs C of dimension 
i is at most 

£ r (=)0PT-< r (=)£l|AVn C '||<r(=)0PT. 

Summing all costs, we obtain the following recursive equation for r(n) 

r(n)OPT < dc d ~ x ■ OPT + d ■ r (-) OPT + r patch (n)OPT . 



Hence r(n) = 0(n max ^ e ' Sc ■•). Choosing c > d l ' £ , as in Step 1, yields 0(n £ ) proving the approxi- 
mation ratio claimed in Theorem [H 

Lemma 10. The cost of patching all relevant cuboids, r patc h( re ), is 0(n £ )OPT. 



Proof. First consider up-patching. Lemma 11 (below) shows the existence of a near optimal network 
that up-patches all relevant cuboids. Lemma [12] shows that by reducing the patching problem to 
gD-RSA, we can find such a network of cost 0(p)OPT, where p is the approximation factor of 
dD-RSA. 

Analogously to the 3D-case there is a approximation-preserving reduction from dD-RSA to DST 



(see Section 2.3), which implies that ciD-RSA is approximable within a factor 0(n £ ) for any e > 0. 

Hence the same approximation factor can be achieved for ciD-RSA by choosing e sufficiently small. 

The lemma follows as the analysis holds analogously for down-patching. □ 

Lemma 11. There exists an M-network of total cost at most (c + l) d OPT that up-patches all 
relevant cuboids. 

Proof. Let J C {1, . . . , d\ be a set of dimensions. For every i £ I we choose a separating plane x % a. 
where ji G {1, . . . , c}. Let J be the set of these separating planes and let C(J) be the intersection 
of C with all separating planes in J. We call C(J) a C-face. There are most (c+ 1) such C-faces. 
Project iV pt onto each C-face. Let N' be the union of all these projections. Clearly, the cost of N' 
is at most (c + l) d OPT. 

We claim that N' up-patches all relevant cuboids. To this end, let (t, t') be a relevant terminal 
pair such that t G Cj lt ...j d , t' G Cj',...,j' an d ji < j[ for all i G {1, . . . , d}. 

We claim that there is an M-path 7Tj from t to Pj 1 ,...j d in N' . To see this, traverse the M-path 
ttu' in iVopt connecting t and t' , starting from t. We assume w.l.o.g. that the separating planes 
that bound cuboid Cj lt __,j d are entered by 7r t t' in the order (xj 1+1 , . . . ,x d , +1 ). The desired path it t 
starts at t and follows it until the separating plane x l - +1 is entered. From this point on we follow 
the projection ntt'( x j -,-fi) of 7r«/ onto C-face C(x 1 - +1 ). If 7r M '(x^ i+1 ) enters £? 2+1 we follow the 
projection vr i4 /(xj i+1 ,^ 2+1 ) of ntp onto C(xj i+1 ,x^ 2+1 )). We proceed in this fashion until we reach 
the C-face C(xj 1+1 , . . . , Xj +1 ), which is just the corner Pj lt ...,j d . Since N' contains the projection of 
7r tt i onto each C-face, the path m described above is contained in N' . This reasoning remains valid 
if the separating planes that bound Cj 1: ...j d are entered in an arbitrary order, as we projected iV op t 
onto each C-face. □ 
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Lemma 12. Given an efficient algorithm that approximates dD-RSA within a factor of p, we can 
efficiently up-patch all relevant cuboids at cost at most (2(c+ 1)) pOPT. 



Proof. In Lemma 11, we showed the existence of a network N' that up-patches all relevant cuboids 
at low cost. Now consider an arbitrary relevant cuboid Cj x j d . Clearly N' n Cj t j d up-patches 
Cj, ,-,. Hence OPTf 3 „• < \\N' n Cj, ,-.||, where OPT" p • denotes the cost of a minimum 

31,— ,3d 31,— ,3d — " 31,—,3d\\1 31,— ,3d 

up-patching of Cj lt .„j d . The problem of optimally up-patching Cj lj .„j d is just an instance Ij 1 ,...j d of 
dD-RSA, in which all terminals in Cj x j d have to be connected by an M-path to Pj lt ...,j d - Applying 
the factor-p approximation algorithm for <iD-RSA to each instance Ij lt ...,j d with Cj ± 4 d relevant, 
we patch at total cost at most 



P 



E °* T Z-,3d £ p E Wnc h ,... Jd \\ < 2 d P \\N'\ 



Cj, : ... : j, relevant Cj. t .., t j, relevant 



The last inequality follows from the fact that each segment of N' occurs in at most 2 d cuboids. The 
lemma follows since \\N'\\ < (c + l) d OPT. D 

Running time. Let T{n) denote the running time of the algorithm for n terminals. The running 
time is dominated by patching and the recursive treatment of slabs. Using the DST algorithm of 
Charikar et al. [3], patching cuboid C 1 - requires time (n!) ( '% where n*- is the number of terminals 
in Cj. As each cuboid is patched at most twice and there are c d cuboids, patching requires total time 
0(c rf )re°' 1 ' e -* = n °( 1 / e ). The algorithm is applied recursively to dc slabs. This yields the recurrence 
T{n) = dcT(n/c) + n ^ < e \ which leads to the claimed running time. 
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